Comment on "Scaling of the linear response in simple aging systems without disorder 
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We have repeated the simulations of Henkel, Paessens and Pleimling (HPP) [Phys.Rev.E 69, 
056109 (2004)] for the field-cooled susceptibility XFc(t) — Xo ~ t~ A in the quench of ferromagnetic 
systems to and below Tc- We show that, contrary to the statement made by HPP, the exponent A 
coincides with the exponent a of the linear response function R(t,s) ~ s~ ( - 1+a * > fu(t/s). We point 
out what are the assumptions in the argument of HPP that lead them to the conclusion A < a. 
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In a recent paper Henkel, Paessens and Pleimling 
(HPP) have addressed the question of the relationship 
between the scalings of the linear response function and 
the zero-field-cooled (ZFC) susceptibility in ferromag- 
netic spin systems undergoing aging after a quench to 
or below Tc- This problem (limited to the quenches be- 
low Tc) had beenpreviously analysed by us in a series 
of papers E3 0] with the following conclusions: 

• the linear response function 



R(t, s) = 



8h(s) 



t > s 



h=0 



scales as 



with 



d — dj 



(1) 



(2) 



(3) 



where z is the dynamical exponent entering the 
growth law L(i) ~ t x / z of the average defect dis- 
tance, dh is the lower critical dimensionality, (n — 
1, du = 3) and (n = 2, du = 4) for scalar or vector 
order parameter, respectively. 

• the ZFC susceptibility xzfc 
scales as 



X.ZFC'(t,s) 



—A 



/*(*/*) 



with 



.4 



a for d < djj 

n/z with log-corrections for 

n/ z for d> djj. 



(4) 



d = d v (5) 



The cases considered by HPP in Q correspond to d < du, 
where the above result gives A = a 6J. Instead, HPP 
reach the different conclusion A < a, thereafter stating 
that A is a new exponent unrelated to a and to aging 
behavior. The purpose of the present Comment is to 
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FIG. 1: Data collapse of s 1 ^ i [xzFc(t, s) — xo] versus x 



show that the data for the same quench considered in 
the HPP paper, if properly interpreted, are in agreement 
with our finding A = a, where a is given by Eq @. 

Following HPP, next to the ZFC susceptibility we in- 
troduce the field cooled (FC) susceptibility XFcif) = 
f dsR(t, s) and the thermoremanent (TRM) suscepti- 
bility pTRM{t, s) — J Q ds' R(t, s'). Obviously, the three 
integrated response functions satisfy the sum rule 



XFc(t) = XZFc(t, s) + f)TRM(t, s). 



(6) 



The HHP argument is built on the behavior of the FC 
susceptibility. From numerical computations they find 



XFc(t) = X0 + Kt 



(7) 



where k is some constant. Making the distinction be- 
tween systems of class S, with a finite equilibrium corre- 
lation length £, and systems of class L, with £ = 00, for 
Xo they make the statement 



Xo = 



(1 




m 2 eq )/T for systems of class S 
for systems of class L 



(8) 



where m eq is the equilibrium magnetization at the tem- 
perature T. Then, assuming that the TRM susceptibility 
obeys the asymptotic behavior 



PTRM(t,s) = s a fu{t/s) 



(9) 




FIG. 2: XzFc{t, s) — Xo plotted against t. 

from Eq. JJjJl follows 

XzFc(t, s)= X o + K>r A - s- a f M (t/s) 
where for the exponent a they take 



1/z for class S 

(d — 2 + rj)/z for class L. 



(10) 



(11) 



Let us now look separately to the different cases. 

1. Class S: 2d Ising model quenched below Tc- 

As an example of class S, HPP consider the d = 2 
Ising model with Glauber dynamics quenched be- 
low Tc- Measuring the FC susceptibilty they find 
that Eq. J2J) holds with xo given by the top line of 
Eq. (JBl and with A = 1/4. Then, assuming that a 
is given by Eq. (|11J1 . that is by a = 1/z = 1/2, they 
make the statement A < a. According to them, A 
does not have any relationship to aging and it is 
due to the roughness of the interfaces, while a is a 
subleading exponent. 

Repeating the simulation of HPP with the same 
quench temperature T = 1.5, we have repro- 
duced their data for the FC susceptibility, However, 
rather than making an assumption on the value of 
a, we have looked for an unbiased comparison of A 
and a, born out of the same set of data. This can 
be accomplished by analysing the ZFC susceptibil- 
ity according to Eq. IjlOfl. In Fig. ^ we have plotted 
s 1 / 4 '[xzFc(t, s )— Xo] as a function of a: = t/s for dif- 
ferent values of s ranging from 1000 to 3000. The 
excellent data collapse in Fig. ^ has two possible 
origins: cither A < a and in the range of s con- 
sidered the third term in the r.h.s. of Eq. IjlOII is 
negligible, or a — A. In order to discriminate be- 
tween these two possibilities it is enough to replot 
the same set of data as a function of t for fixed s. If 
the first alternative is the right one the data should 
collapse also in this plot, while in the second case 
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FIG. 3: XzFc(t, s) — Xo plotted against s for different fixed 
values of x. Straight lines are the best fits (with exponents 
0.24,0.25,0.25,0.27,0.27 from top to bottom). 



there could be no collapse, due to the existence of 
the s dependence. Fig.|21shows that indeed the lat- 
ter one is the case, that the s dependence is a large 
effect and that, therefore, a = A = 1/4. 

Although the pair of Figures ^ and [21 is certainly 
enough to settle the issue, in order to help visualize 
the result we have also produced (Fig.|3) the log-log 
plot of xzFc(t, s) — xo as a function of s, for fixed 
x. The data in Fig. [3] display an excellent fit with 
a single power law, with a slope very close to —1/4 
for every value of x, which makes clear at glance 
the conclusion a = A = 1/4 reached above and in 
agreement with Eqs. J2| and JSJ). 

In summary, this means i) that a is not subleading 
as believed by HPP and ii) that roughening of the 
interfaces, rather than being unrelated to aging, is 
precisely the mechanism that renders the exponent 
a smaller (for d < djj) than 1/z, as explained in 
ref. . What goes wrong in the HPP interpretation 
of the data in [l| is the assumption that, for systems 
of class S, a is given by the top line of Eq. I|llfl 7]. 

2. Class L: 3d spherical model quenched to and below 

As an example of class L, HPP consider the spher- 
ical model quenched to and below Tc- They com- 
pute numerically XFc(t) in both cases, with d — 3, 
finding that it saturates to a constant for large t. 
However, they do not identify this constant with 
Xo, since they make the statement (|SJ| that xo = 
for systems of class L. This is supported by an hand 
waving argument according to which xo ought to 
vanish, since the correlated clusters for systems of 
class L should have no "inside". Although vague, 
this assumption is crucial because it is the starting 
point of the chain of implications used by HPP: the 
saturation of XFc(t) to a constant value and xo = 



imply A = 0, which in turn implies A < a, since 
in the spherical model a = (d — 2)/2, both in the 
quench to and below Tq. 

Before going further, it is necessary to clarify the 
physical meaning of the constant xo appearing in 
Eq. Q. This can be readily understood recalling 
that 

lim X Fc(t) = Xc, g = (1 - m 2 )/T (12) 

where Xeq is the static susceptibility. Since the 
second equality in the above equation is nothing 
but the static fluctuation-dissipation theorem, it 
holds for all systems, be they of class S or class 
L, quenched below Tc or to Tc- Therefore, from 
Eqs. and ljT2)l. we can identify xo = Xeq aim 
Eq. JSJ must be replaced by 

Xo = (l-ml q )/T (13) 

both for systems of class S and class L, which 
clearly reduces to xo = for the quenches to 

Tc- Therefore, Xo = assumed by HPP is excluded 
in all cases. 

Furthermore, in the case of the spherical model 
there is no room for assumptions, since the model is 
exactly soluble. Using the formulas that HPP give 
in 0, it is not difficult to derive analytically the 
large t behavior of the FC susceptibility obtaining 

XFc(t) = (l-m 2 eq )/T + Kt-^/ 2 (14) 

with m 2 eq = 1 - T/T c , K = [r(l - d/2)r(l + 
w/2)]/r(2 — ^r-) where T is the gamma function, 
lu = d/2 — 2 in the quench to Tq and lj = —d/2 in 
the quench below Tc- Therefore, comparingwith 
Eq. 0, we have xo = Xeq = (1 - rn 2 eq )/T [U, as 
expected from Eq. I|13fl . This implies A — a — 
(d - 2)/2 for T < T c , in agreement with Eqs. Q 
and ©. 

3. Class L: 2d Ising model quenched to Tq- 

As an additional instance of a system of class L, 
HPP consider the 2d Ising model quenched to Tc- 
Again, they find that the FC susceptibility sat- 
urates to a constant. By the same reasoning as 
in the previous case, from the assumption that xo 
ought to vanishes they make to descend A = < 
a = {d — 2 + rj)/z c = 0.115, where we have used 
z c = 2.167 and the exact result n — 1/4. 

Although the argument of Eq. (|12fl ought to suf- 
fice, we have repeated their simulations and we have 
plotted (Fig. 2} log(xo — XFc(t)) against logf, with 
Xo = l/7c. The figure shows a very clean power 
law decay with A = 0.115 ± 0.005, which compares 
very well with the value of a given above and yields, 
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FIG. 4: The field cooled susceptibility XFc{~t) — Xo vs t for 
the 2D Ising model at T = T c ~ 2.269. The straight line is 
the expected behavior XFo(t) — Xo ~ t v ^ Zc with r\ — 1/4 and 
z c = 2.167. 

again, A = a. Furthermore, the observation of the 
decay with the correct value of the exponent im- 
plies that also the subtraction by xo = 1 /Tc is the 
correct one. 

In summary, the unbiased analysis of the data for the 
FC susceptibility in all cases considered by HPP, that is 
in systems of class S and class L with d < du, yields 
A = a in agreement with Eqs. and JSJ. The biases 
in the HPP analysis, which lead to the wrong conclusion 
A < a, are in the two assumptions i) a = 1/z for systems 
of class S and ii) Xo = for systems of class L. 
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